Appendix A 

Derivation of mutual inductance 
between two dipole coils 

Hero. « derive the mutual inductee between two mastic dipole coil* 

where , A 

a = unit vector in direction of dipole .A 
* = magnitude of effective area of dipole A 
% = unit vector in direction of dipole B 

B = magnitude of effective area of dipole O g 
r = unit vector in direction from center of dmo e A to cent P^ ^ 

R = magnitude of distance from square coil. Each coil Ls shrunk 

Each dipole is replied by a smaU but fim^ ^ obtain the dipole bim t. 

to an infinitesimal size (while ^^^^Z L «** properties are described by 
In the dipole limit, the shape of ^ " * trical area and number of turns. 
the coil's effective area, the P/^^^^ts is given by Neuman's formula: 
The mutual inductance, L m , between 

-if. i / (A - 2) 

Um ~ 47T /circuit.* Ja-rmitD r AB 

distance between point on circuit A and point on circuitB 
Define some quantities: 

■4 = = effective area of dipole A ^ ^ 
fx and v are two unit vectors such that a = u X v 
# = 23b = effective area of dipole B ^ 
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Use the S, 0, a cartesian coordinate system, where dipole A is at the origin and dipole B is at 
~R = Rr = 7xu + vv + aa. 

The four corners of coil A are: 

^2 = 1^- 

Not* th!t U the'ratio of the effective area of dipole A to the effective area of coil A is 
A point on coil A is given by: 

where p. has nothing to do with 
The four corners_of coil B are: 
Bl=k 2 p + %q + R 

B~ 2 = %p-\q + R_ 

Not" that the^ratifof the effective area of dipole B to the effective area of coil B is £. 
A point on coil B is_given by: 

(p-.<?) = PP + <IQ + R ., „ . 

The distance from a point, on coil A to a point on coil B is: 

tab = \sA-pp-n<i I . .... 

where the vector from center of coil B to the point on coil A is: 

B^dtVil^ around_circuit B in equation A.2_by explicitly going around circuit B from 
point Bl through points Bj, B 2? B 3) and back to pomt B 4 : 

u B r + ff Pdb Qdb P db , ^ -\dsA 

' ' (A.3) 

Define g to be the inverse of the distance from a point on coil A to a point on coil B: 

1 1 (A.4) 

Calculate some directional partial derivatives along the p and q axes, remembering that sa is 
independent of j> and q: 

\dqV = * d r 1 1 -lim— f - r ^ i-r l (A.5) 

[SET = I fcd^U--* = {w^m ~ wrwrm) (A7) 
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Set e = h, and substitute into equation A. 3": 

(A.9) 



Perform the integration over b: 

Expand the integral around circuited in equation A 10 by explicitly going around circuit A from 
ttt nnints X". TG, A 3; and back to point Af. 



where: 



t~ = ! ^S3,( / ' + /! + i3 + /4) 


(A.11) 


: ^ / ( ? [|] ,-,.0 " ? til „» ' 8 * 


(A.12) 


/ H%]^=o~ ? [^L,=o) M =|* r,dl/ 


(A.13) 


_ 2 
"2 


(A.14) 




(A.15) 


the four dotproducts: 




/x^U (/ . 7 . S _ / 5 . G - V p • 5 + 9 • v) 
4tt h-o v 


(A.16) 



-2. f (r*r* 

-J 



where: 

j * i it"* ^ i i — — t i u-/a (A-'^-'O 
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(A.18) 


^-tfj \Ud 7;=9= „ L&Z 

h 


Jp= f/ =0/ 


(A.19) 






(A.20) 


Apply the definition of a partial derivative and the fact that h is small to equations 
A.19, and A.20: 

i Ha 2 9 r , a5 l 


A.17, A.18, 
(A.21) 


2 


f a 2 . V i= " =0 
L^Jp=,=o 


(A.22) 


2 r .-,0 n/i=^> 
2 


[ a 8 , 

|/W<zJ p=<?=0 


(A.23) 


'-•-hi [o^ep\ p .^ 


r gig lt*= u=0 


(A.24) 



Recall that g is: 

! 1 1 _ ^ ^ (A.25) 

^ _ = _ (A-26) 

_ [_^_r =v=o = - (a - 27) 

, f^r- 0 __[j^n =0 



(A.28) 
(A.29) 
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Substitute into equation A. 16: 

u. AB ( d 2 g „ ~ &9 ~ ~ , d 2 9 s . d _^_S-v\ (A.30) 

A <m and its partial derivatives in the 2, v, 5 cartesian coordinate 
Calculate g in equation A 30 and its paraaj ^ 
system, where dipole * is at the origin and dipole B » at A 

J^L = 3g 5 (p ■ un + P ■ vv + P • a«) « ~ -9 3 P " G 
= 3 5 5 (p • uu + P ■ vv + p • So) « - 3 3 P • « 

- 3o 5 (o - Gtt + 0 • vv + q- aa) » ~ 0 3 9 ' 2 
Out), oy v ^ 

= 3j 5 (?-®»* + ?- 8to + ?" oa ' , '~ flrfl " 1 ' 
Substitute into equation A.30: 

Lm_ 47rR 3 VH 2 / 

where _ . r\ 

T = 2 ((p - v)(S -S) — (p*«)(g-v)) = ~2(a • b) 

v u + (p • fit* 4- p • vv + P • 2a) q ■ v u 

Insert now cancelling terms at end: m™_(f-i!l(S »la , + (? «)(?• 

Since o = u x v and ft = p x g, we have: 

-s-6 = (p«)(5- fi )-(P* s ] (if ' 5J) , • ~ 
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a)(q v)ua + (a b) R? „ 

Since o = Sxvandft = px 9l w> have: 

^ * (ggKjg) vzd&m «- ! fiBfia "'-^^ «j±<& m 

(a • f)(i> ■ r) = « 

Then: ~ ~ N „, 

p = -(a - f)(b • f) i? 2 + (S • b) -R 
Substitute F and T into equation A.31: 



Lm = ( 3 (2 • r)(b ■ r) - 3(8 • 6) + 2(S - 5)) 

4.7T it* ^ 

Simplify to get result A.1: 

We can write equation A.33 using a dyadic operator: 

Lm 47T R? 



(A.32) 



(A.33) 



(A.34) 



If the vectors a, I, and r are three-element row vectors then the dyadic operator V can be 
written^ a 3x3 matrix, and the dot product, are matnx mulUpbes: 

(A.35) 



where V is given by (I is the 3x3 identity matrix): 

V = 3f*f - I 



(A.36) 
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Appendix B 

Derivation of orthogonal-dipole seed 
calculation 

Given ». -del of mutual inducts L, - «— - 75 M " 1 o^ 0 " °' - 

^.^rrplo ^dO. seed, we cnicniate -n— vames of S and O 
""^n^Tc tf£e apprcedrnation, that the dimensions of each eon are sma.! compared to the 
distance between coils. ffArifll j:„ rt ip transmitter coils and three equal- 

diX in Appendix A, so we start with a 3x3 mutual-inductax.ee matnx: 

Ltr = f^\z(t.R)(R-r)-(tr)) (B.1) 

where 

t = transmitter coil 0J. 5 or 2 
r = receiver coil 0,1. or 2 
/x 0 = 47rl0~ 7 henries J meter 

R = magnitude of distance from dipole t to dipole r 

R = unit vector pointing from dipole t to dipole r 

4 t = magnitude of transmitter dipole effective area independent of t 

4 = magnitude of receiver dipole effective area independent ot r 

't 1 unit vector pointing in direction of dipole t effective area 

f = unit vector pointing in direction ofdipole r effective area 

The transmitter effective area vector A t is: 

Ai = A t t . 

The receiver effective area vector A r is: 

~A~ V = A r r 

Define the cartesian coordinate system: 
X = (1 5 0,0) unit vector in +X direction 
Y = (0,1,0) unit vector hi H-V direction 
2 = (0,0 : 1) unit vector in +Z direction 
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The receiver position vector ft is independent of t and r. and has cartesian coordinates fix, B. . 

and ^ ^ 

ft = ftft = RxX + RyY + , 2 which MC centered on the origin, orthogonal to 

We have three transmitter dipoles, 0, 1 and <s, Tnuicn cu respectively: 
each other, all the same effective area A, and pointing along the A , Y, and Z axes resp, 



Wheat = 0, t = X. 
Whent= 1, t = y. 



Whent = 2, receiver dipoles. to represent the receiver 

Define a rot ated cartesian coordinate system for the recener cup 

orientation: 

Xr = unit vector in rotated +X direction 
Yr = unit vector in rotated +Y direction 

f r = unit vector in rotated direction mtation Quaternion O. Then we calculate 

complement of O): 

(0X) = 0(0,1,0,0)0' 
(O.^) = O(0, 0,1, o)o' 

(0,Zr) = 0(0,0,0,1)0 on the int R> orthogonal to 

We have three receiver dipoles, 0, 1 and l, wnicn a respectively: 
each other, all the same effective area Ar, and pointing along the X rt Y r , 
When t = 0, f = A r . 
When r = 1, f = Yr. 
When r = 2, r = Z r . ■ _ ^ 

Define A*, which is independent of ft, O, and Ltt . 

Now B.l becomes: K v fB 2 \ 

L tr = ^(3(t.ft)(rfi)-(^0) ^ 

gomB^ 

L?r = ( 9 (t • ft) 2 (f • ft) 2 - 6((t • ft)(r • fi)(t • F)) + (t • f) 2 ) ( B3 > 
Then, the sum of the squares of the mutual inductances of one transmitter clipole and the three 
receiver dipoles, is: 2 

L?-|r(8(r.«» + i) 

For transmitter coil 0: ^ - , (B 5) 

L ° = ft« V V 

For transmitter coil 1: 2 K* ( , f Ry, 2 . A (B.6) 

^ = ^ V 3( fi" ) V 

For transmitter coil 2: x2 „ , (B ?) 

1,2 = ftS" V 3t ft J V 
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The s\irn of the squares of all nine mutual inductances is then 

K' 2 

'7? 



lit = (b .8) 



Solve B.8 for R: , 



Solve B.5, B.6. and B.7 for Rx 2 . Ry 2 and J?.z 2 respectively, and use B.8 and B.9 to eliminate 

*-WH)ffi)' 



Now we take the .square roots of these three quantities: 



The signs of Rx, Ry, and fire determined from the signs of products of the mutual inductances. 

Due to the symmetries of the fields, positions +R. and —R, for the same orientation O, give 
the same mutual inductances, and cannot be distinguished by calculation. In the literature, this is 
called the hemisphere ambiguity. 

To resolve this ambiguity, we choose to use the + X hemisphere, where i?* is always positive. 

Form sums of products of mutual inductances L tr : 

S 0 i = LoqLiq + LqiL u + ^02^12 (B.16) 

520 = ^20^00 + L21L0} + ^22-^02 (B.17) 

Recall that the signum or sign function is denned to be: 

sgn(x) = +lfor{x > 0} ? -l/or{> < 0} 

Then the signs of Rx. Ry, and Rz are given by: 

sgn(Rx) = +1 

sgn(R y ) = sgn(Soi) 

$gn(Rz) = sgn{S2o) 
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Then the position R cartesian components are emulated in terms of the mutual inductances 
L * by: Ll = Woo + LoiLo, + UnU* < B - 18) 

r T (B.19) 

L? = LioLia + LnLn + 

Ll = L^L-io + L 2 i^2i + £22 £22 

i (B21) 

[Tig A /6J£\* (B.22) 

, r , . ff 2 L « - ±\ f^V (B.23) 
^ = sgniLooLw + L01 L n + £02^12) ^ ^ ^ 3 j V I&t / 

r r r 1 I f 2 L '° - (B.24) 
= sgniLxLoo + ^21^1 + ^22^02)^ ^j^r 3 J V L£* / 

M H„, and/or ft* ^ ^ E 

t rn^si^ss. - ~ ^ — ^ ^, . ca„» 

a close-to-zero result to become imaginary. ,^ Watpd dioole moments combine the same 

Since the transmitters are colocated dipoles and ^Transmitter dipole triads 

way vectors do, synthetic mutual-induct.^ce ma^s ^ numericaJ instabmty of 

with other orientations than along the X, "r , and Z axus. 

the above solution. , , h ■ numerically unstable), we 

If the calculated position is close to the X, Y, or Z axis calculate 

Calculate unrotated position R" = (K T , /V ^ oy suds™ & 
tances L tr in B.22, B.23, and B.24: 

, r x I^A-I^ (B.26) 
R% = afl n(Loo^io + Lox U 1 + ^02^12) ^ ^ ^ zJ\L^J 

= s S 7x(L2 0 Loo + ^-21^01 + 1-22^02)^ ^2-^- - -J \-£T) 

if of ff» or R a come out imaginary, make that value or those values zero. 
If any of R x , i^, or K s come 01 b - this causeg olJy a sma n 

Note that the sign for a small-value f£, fq,, or K z may dc wi« 5 
position error, which error will be removed in following steps. 
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Calculate a unit vector R u in the direction of R": 



R u = 



(B.28) 



-( 2 



y?(/?S) 2 -K^) 2 + (tt!) 2 ) 

We next calculate a unit vector 5 s in the (* fl n(i?£), sgn(R^), sgn(R%j) direction: 

■gn(Rl) sgn(RX) sgn(R%) \ (B 29) 

Note that 5" is far from the X, Y, and Z axes, and is alwaysm the +X hemisphere. 

Note that our choice of ^ means that the dot product. R? • S» can never be negative. 

Note that R u ■ & analytically can never be greater than unity. Numerical errors can result in 
a value slightly greater than unity, but such errors do no harm here. 

If & ■ 5 s is greater than or equal to some threshold (we use 0.9 to 0.95), then the position 
estimate is numerically stable enough to be accurate, so the position is given by: 

R = 7F (B30) 

and we skip equations B.31 through B.45 below. 

If the dot product is less than the threshold, then the position estimate is mimencally unreliable, 
and we must improve the numerical stability. We do this by recalculating the position in a rotated 
coordinate system, and then derogating the position back to the original coordinate system. 

We choose to rotate JZ» into 5". Note that the rotated position vector will be far from the X, 
Y and Z axes, even if the recalculation changes the position a little. 

? We could use any of various vector-rotation formalisms. We choose to use normalized rotation 
quaternions for their numerical stability and lack of singularities. 

The angle between R" and S u is: 

9 = arccos(R» - 5") _ _ 

The normalized rotation quaternion to rotate into S u is: 



U = ^cos (?) 



R* x S u 



cos I 



(!) 

and 



(B.31) 



(0, S^)=U (0, H") W (B-32) 

Equation B.31 is numerically unstable when R» and point in close toopposite directions, i.e. 
when | approaches ninetv degrees. This will never happen here, because R» ■ is never negative 

We treat the 3x3 matrix of mutual inductances, Ltr, in B.2 as a trio of column (transposed row) 
vectors, one for each receiver coil, where the superscript t indicates the transpose operation: 

L lr = (LU'i.4) ( R33 > 

To rotate the transmitter coordinates, we apply the rotation U to each of the three row vectors 

in B.33: . 

(0,I^)=U(O.L 0 )W (B-34) 
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(0,15)- W(0. ia)W (B 36) 
The rotated mutual-inductance matrix is: 

/ rr t jrt TTt\ (B.37) 

, , m - (W R r R r ) by substituting the rotated mutual inductances 

Calculate rotated position R r = \K x ,K y ^z) & 

L lr in B.18 through B.24: ^ = ^ + ^ + ^ (B.38) 

==^0^0+^21 + ^*2 C • ' 

LL = ^ + L ? + LI ^ <B*> 

1/ Lg 2 A (B.42) 

[7^nT\ f 6JCf2 V (B.43) 
Rl = 8jn(LooLio + £oi L " + >| ^ 2 J?^ Z)\L 2 tat ) 



\ ( 0 L% _ A (B.44) 

position is far from the axes. 

Derotate the position: = ^ ^ ^ (B.45) 

When r = 0 ? f = * r 
When r = 1> f = 

When r = 0, r° = X. 
When r = 1. r° = 

Thf orientation ns . nom.lfceo rotation ,uatemion O. 
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O (o.y)o' ^ 

We can afco the oriental - - ort,_l — matrix O. « 

the vectors are row vectors: _ ^ (B.49) 

(B.50) 



Y r = YO 

— ~„ (B.51) 
Z r = ZO 



(indicate., the transpose of a vector or matrix: ^ ^ 

*-(*«.«*.«')(*. W (B ' 53> 

matrix at position J?, for zero orientation: 

L O r = ^(3 ( T.R)0-.R)-(frO)) ( B ' 54 > 

L? by the the receiver rotation: (B.55) 
Solve B.55 for O: _ ^ £0 -i (B-56) 

O U - or— ~*^Jg~ 

best we can do. is to convert lo . be** eSe^n t^X^terUs' , by Berthold 
^^SStX^ SrSSS America", ^ 4. Aprii, .OS, p*. 
629£f , herein incorporated by reference. 
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Appendix C 



Mutual inductance of circuits of 
straight-line segments 



The mutual inductance, L m , between two closed circuits is given by Neuman's formula: 

TcircuilA JcircuitB TAB 



L m = ^i I mm (CI) 

4/T J circuit A /circuit 



(C.2) 



where 
and 

u„ = 4tt 1(T 7 henries /meter 4. 
Here, we consider the case where both circuits are composed of straaght-lme segments. 
Circuit A is composed of segments A\ • ■ - A n • 
Circuit S is composed of segments By • ■ f? m . 
Then C.l becomes: 

4t / Bf*» rt Aj 1 -"" *VlB 
The integral over A k can be evaluated analytically: 

where 

s start vcc tor from origin to start of segment Ak 
s{ inish = vector from origin to finish of segment A h 
s ^ nrt = vector from origin to start of segment Bi 
s[ inUh = vector from origin to finish of segment B\ 

— r f itusfx start 

Ti = s{ iilish - sf Qrt 



4 = ^ ■ Sk 
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sf = si - si 



T = (sf tart - s s k art ) 

s s g = +i if o), -i£pc-^< o) 

_ = j f (Sfl = +1)> sf t*rt if (Sp = -1) 

sj; = *f tar * - 37 + sib 

* k Z Si = cosine of angle between S6 and Sfe 
1 7*FI . w-li mimer icallv by well-known numerical integration methods. 

Equation C.3 is easily evaluated mm ™™ l *j£J™ . opposite directions along segment A k . 

C 3 - r ■ i • ei..Mi»v nrnirs when 37 and approach being antiparallel. 

The second region of numerical instability ^ xa ^ h ^/' £ become exact i y parallel, 

in this region, * • « < 0, and 1 + c t approach « « ™£ ^ denominator in the integral 
1 + ct = 0 exactly. In this region, we choose s g - -1 to avoia a 

" Td^eUion C3, start with the term in equation C.l for one pair of current segments: 

,flf*** M?'"* h m-dSB (C.4) 

Current segments * and * are both straight lines so it is natural to describe them in linear 
terms of dimensionless variables a and & respectively, where 
0 <a < 1 

0 - h ~ 1 ! „i;™1W first ^ we will keep dependendencies on a explicit. Depen- 

7 — r finish 

When % = +1, then < A (fl) = ft -*fc° 

When s 9 = -1, then M (a) = -jf* + »*» , . 
A point on segment B, is a linear function of b: 

J^) = sf art + sib 
Then: 

(tsB =~sidb 

*s-j.d tt= -s„7b ■ si da db 

r AB = y/(*b + 5 9^°) " S + 
Now eciuation C4 becomes: 

yl r(i-s 9 )/2 da ^ (C.5) 
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Expand the denominator in equation C.5, separate out the coefficient of a* ; and remove the s 9 
dependence from the limits of integration: 

= ^7 f l f 1 da db 

W y/4 ^ ° ^/f| + H£alpi a + a 2 

Define two quantities which do not depend on a: 

•Note that 0 < 8 < 2. 
Define a new variable x: 
x = a + a 

Now equation C.6 becomes: 

We have for /? > 0: 

/ / 3 = log(x + + 0 (C - 8) 

Substitute equation C.8 into equation C.7: 



'« = io log( oT+ToTO? 



Substitute the definitions of a and 0 into equation C.9: 



, /I? Jo // /*T . \2 _l fin — r?l 



Remove a common factor of from numerator and denominator: 



Sk ' si r w / v * V Y )db (*■/. 

77 Jo + v^^To^i) 



Simplify: rr /„. 7-755T^i 



Jt ,^/; toe r + ^;V +i -^_ ) , 
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Substitute for c* in the last radical in the numerator: 



^2 Jo ° K S 9 C t +l 



T M = / log( , Y , ; ' ) dft (C13) 

Factor the numerator inside the last radical in the numerator: 



lH= S ±4 log( V » ~* — r ) db (C.14) 

2 Jo 5 p Cf + 1 



s i_ 

Substitute equation C.1.4 into equation CA, and that into equation C.2 



This is our result, equation C.3. 
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